Forecast Constraints on Inflation from Combined CMB and Gravitational Wave 

Direct Detection Experiments 
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We study how direct detection of the inflationary gravitational wave background constrains infla- 
tionary parameters and compliments CMB polarization measurements. The error ellipsoids calcu- 
lated using the Fisher information matrix approach with Planck and the direct detection experiment, 
Big Bang Observer (BBO), show different directions of parameter degeneracy, and the degeneracy is 
broken when they are combined. For a slow-roll parameterization, we show that BBO could signif- 
icantly improve the constraints on the tensor-to-scalar ratio compared with Planck alone. We also 
look at a quadratic and a natural inflation model. In both cases, if the temperature of reheating is 
also treated as a free parameter, then the addition of BBO can significantly improve the error bars. 
In the case of natural inflation, we find that the addition of BBO could even partially improve the 
error bars of a cosmic variance-limited CMB experiment. 
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I. INTRODUCTION 

Inflation [l|-[3j], which is widely believed to have taken 
place in the very early universe, generically predicts the 
existence of tensor mode perturbations originating from 
quantum fluctuation ■ Since inflation predicts an al- 
most scale-invariant spectrum, the tensor mode pertur- 
bations are considered to exist as a gravitational wave 
background over a wide range of frequencies. Because of 
their weakness of interaction with matter and their lin- 
earity, they are expected to remain uncontaminated even 
at higher frequencies, whereas the scalar mode perturba- 
tions become nonlinear and it might seem difficult to re- 
cover information about its spectrum on smaller scales. 
Therefore the gravitational wave background possesses 
the potential to enable us to reconstruct the primordial 
spectrum over a wider range of scales. This helps us learn 
more about the inflation potential since the amplitude of 
the spectrum at each scale corresponds to the height of 
the potential when the mode exited the horizon during 
inflation. 

One of the ways to detect the inflationary gravitational 
wave background is to measure the i?-mode polarization 
of the cosmic microwave background (CMB) induced by 
primordial gravitational waves @, [1] • The ongoing satel- 
lite mission, Planck |9j, is expected to detect the signa- 
ture of the large-scale tensor perturbations if they are 
sufficiently large in amplitude. The other way of probing 
the gravitational wave background is to detect it directly 
with space-based laser interferometers such as Big Bang 
Observer (BBO) and DECi-hertz Interferometer Gravi- 
tational wave Observatory (DECIGO) (MCLTI- Clearly, 
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both approaches provide us with new information about 
the early universe, and help in determining inflationary 
parameters more accurately if they detect the gravita- 
tional wave background. It is notable that while CMB B- 
mode experiments measure the gravitational wave back- 
ground at the present horizon scale (~ 10 -18 Hz), direct 
detection experiments measure it at much smaller scales 
which correspond to the detector size (~ 0.1 — 1Hz). This 
means that these two different types of experiments could 
provide independent information. 

In this paper, our goal is to predict how the direct de- 
tection of the gravitational wave background can com- 
plement CMB observations. The complementarity of 
the two observations is discussed within the slow-roll 
paradigm in Refs. [l2l-[p3]. While they only discuss the 
detectability of the gravitational wave background in di- 
rect detection experiments by connecting the amplitude 
of the gravitational wave background at direct detection 
scales with model parameters which are allowed by the 
current CMB constraints, we perform the Fisher matrix 
calculation and forecast errors on the parameters attain- 
able in the direct detection experiments. This enables 
us to discuss how well inflation parameters may be de- 
termined by direct detection of the gravitational wave 
background and how degeneracy in parameter space may 
be broken when the two constraints from CMB observa- 
tion and direct detection are combined. 

Here, we consider the combination of the upcoming 
CMB satellite experiment Planck and the future satel- 
lite gravitational wave detector BBO. As long as we 
make the usual assumption that the equation of state is 
> — 1 during inflation, the spectrum of the gravitational 
wave background cannot be blue-tilted and sensitivities 
of ground-based experiments and the preceding space 
mission LISA are not enough to detect the inflationary 
gravitational waves. Therefore, it would be a long road 
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to seek the inflationary gravitational waves, while CMB 
experiments are expected to detect -B-mode signals in the 
not so distant future. One may think constraints on in- 
flationary parameters may be improved by post-Planck 
CMB experiments before direct detection is achieved by 
BBO. For this reason, we also look at the case combined 
with a cosmic variance-limited experiment, which would 
be similar to CMB satellite missions planned for a few 
decades like BBO, instead of Planck. 



Applying the Fisher matrix method to the specifica- 
tion of Planck and BBO, we evaluate errors on infla- 
tionary parameters expected in these two future experi- 
ments. To calculate the Fisher matrices, we need to ob- 
tain differentiation of the spectrum with respect to the 
model parameters assuming a fiducial model. The CMB 
spectra are computed using the CAMB code [15j . For 
the gravitational wave spectra at direct detection scales, 
we present two different approaches: an analytical way 
with the slow-roll approximation and a numerical way. 
We first apply the Fisher matrix method to the slow- 
roll paradigm, which enables us to express the spectrum 
in terms of the parameters defined at the CMB scale. 
This slow-roll framework is a simple and well-established 
way to connect the amplitude of the spectrum at differ- 
ent scales. However, since this method uses the Taylor 
expansion to extend the spectrum from the CMB scale, 
the approximated spectrum may deviate slightly from 
the true value at the direct detection scale [Tg] . For 
this reason, we also present the case in which numerical 
calculations are used to obtain the precise amplitude of 
the spectrum. This method, however, requires us to as- 
sume an inflation model. Therefore, both of the methods 
have their advantages and disadvantages; the slow-roll 
paradigm can be applied as a more general model, which 
does not require us to assume an inflation model like the 
numerical method, while the numerical calculation en- 
ables us to make a more accurate prediction. 



The outline of this paper is as follows: In Sec. UH 
we briefly introduce the Fisher matrix methods for both 
CMB measurements and direct detection of the gravita- 
tional wave background. In Sec. IIII1 the Fisher method 
is applied to the slow-roll paradigm. We show how direct 
detection reduces errors obtained from CMB in the ns — r 
plane by connecting ns and r with parameters of the 
gravitational wave spectrum assuming slow-roll inflation. 
In Sec. IIV1 the investigation with numerically calculated 
spectra is given for several inflation models which predict 
a large enough amplitude of the gravitational waves to be 
detected by BBO. We also allow the temperature of re- 
heating to be a free parameter. We show to what extent 
potential parameters can be determined when the con- 
straints from CMB and direct detection are combined. 
Conclusions are given in Sec. |Vl 



TABLE I: Planck instrument characteristics 



Center frequency (GHz" 


) 70 


100 


143 


217 


9 (FWHM arcmin) 


14 


10 


7.1 


5.0 


<7T (flK) 


12.8 


6.8 


6.0 


13.1 


CJE (^K) 


18.2 


10.9 


11.4 


26.7 


a B (fJ.K) 


18.2 


10.9 


11.4 


26.7 



II. FISHER MATRIX METHODS 

The Fisher information matrix is commonly used in 
cosmology to predict how well parameters can be deter- 
mined in future planned experiments, and is defined as 
the second derivative of the log likelihood function C with 
respect to the model parameters Pi, 

According to the Cramer-Rao inequality, the Fisher ma- 
trix gives a lower bound for variances of the parameter 
estimates. This enables us to estimate expected errors 
on model parameters for a given experiment. The like- 
lihood function is constructed from the specifications of 
each experiment. 

A. CMB B-mode polarization 

The CMB Fisher matrix is given by (see for example 

^-EEf^W^f)^- (2) 

i x,x> ° Vl op i 

where X and X' are summed over the the CMB temper- 
ature, i?-modc, and B-mode of the CMB polarization. 
The covariance matrix can be obtained from Zaldarriaga 
et. al. [l7j and depends on the temperature noise per 
pixel (or), the polarization noise per pixel (<7e and erg), 
the pixel area in radians squared (0 2 = 4?r/iVpi x ), and the 
beam window function which we approximate as Gaus- 
sian (Be. ~ exp(— £(£+ l)of). The values we use are taken 
from the Planck blue book 1 and are listed in Table Q] 
(note that 9 needs to be converted to radians). We use 
<7b = 9/ \/81n2 and combine the different frequency bands 
as specified in Bond et. al. [ill]. We take the range in I to 
be 2 to 2000. At higher £, secondary sources of temper- 
ature and polarization will likely prohibit the extraction 
of cosmological information from the primary CMB. We 
assume that the foregrounds can be removed by using 
templates and multifrcquency information; see, for ex- 
ample, Efstathiou et al. Efstathiou et.al. We also 



http: / / www.rssd.esa.int/SA /PLANCK/docs/Bluebook-ESA- 
SCI(2005) l_V2.pdf 



3 



consider the case of a cosmic variance-limited (in both 
temperature and polarization) CMB experiment (CV) for 
£ < 2000 and without delensing. This would be similar 
to optimistic foreground removal with epic-2m (20| . 

Throughout this paper, we assume a flat A cold dark 
matter Universe and we use the WMAP5 maximum like- 
lihood values for the nonprimordial power spectrum pa- 
rameters [H|: (baryon density) fl b h 2 = 0.0227, (CDM 
density) Cl c h 2 = 0.108, (amplitude of curvature perturba- 
tions) A 2 lm = 2.41 x 10~ 9 , (reionization optical depth) 
t = 0.089, and (the Hubble parameter) h = 0.724. 



where the wavenumber relates to the frequency as k = 
27r/. The critical density of the Universe is defined as 
p c = 3H 2 /8ttG, where H is defined by the scale fac- 
tor a(t) as H = a/a, and the energy density of the 
gravitational waves pew is given by the 00 component 
of the stress-energy tensor. Let us consider the tensor 
perturbation in a Friedmann-Robertson- Walker metric, 
ds 2 = —dt 2 + a 2 (t)(Sij + hij)dx % dxK It is convenient to 
expand hij into its Fourier components, 



lij (t , x) 



E 

A=+,x 



d 3 k 
(27T) 3 / 2 ' 



Ak)hi(t)e 



ik-x 



(7) 



B. Direct detection of gravitational waves 

In a direct detection experiment, cross correlation 
analysis is a powerful method to detect a weak stochastic 
gravitational wave background, such as inflationary grav- 
itational waves [22M25I ] . The Fisher information matrix 
for the cross correlation analysis is given as 26] 



*E 



3Hl 
10n 2 



df 



2T< 



obs 



hu(f)\ 2 d Pt n GW (f)d P] n GW (f) 

/ 6 5/(/)5j(/) 



(3) 



where / is the frequency of the gravitational waves, H 
is the Hubble constant, T b s is observation time. Here, 
we consider TDI (Time-Delay Interferometry) channel 
output (7 = A, E, T) which would be adopted in the 
BBO project. In this case, the noise transfer functions 
Sa,e,t(J) are given as [27H2I 



where the polarization tensors e^' x satisfy symmetric 
and transvcrsc-tracclcss conditions and are normalized 



as e ij( e ij)* = 2<5 AV . Then p GW is given as 



1 



((5 t ^) 2 + (-V^-) 2 > 



64ttG 
1 

32ttG J (2tt) 3 a 



d 3 k k 2 

3 72 



2Ei^i 2 , 



which yields 



CI 



GW 



1 / k 

12 \aH 



Ei fc 



A 1 2 
kl ■ 



(8) 



(9) 



One may use the tensor power spectrum A^(fc) instead 
of f^GW, which is defined as 

dihijhV) 



din k 



E^ 



A 1 2 
kl ■ 



(10) 



S A (f) = S E (f) = 8sin 2 (/72)[(2 + cos/)^(/) 

+2(3 + 2 cos / + cos(2/))^(/)], (4) 

StU) = 2[l + 2cos/] 2 K(/) +4sin 2 (//2)^(/)], (5) 

where / = 2irLf and L is the arm length of the detector 
which is assumed the same for each arm. In the case of 
the standard BBO detector, the arm length is L = 5.0 x 
10 4 km and noise functions are S*£ = 2.0 x 10" 34 /L 2 Hz _1 , 
S% = 9.0 x 10- 34 /(2tt/) 4 /(2 j L) 2 Hz- 1 . 2 The overlap 
reduction function "fu(f) is calculated with information 
about relative locations and orientations of detectors (30l — 

n. 

The intensity of a stochastic gravitational wave back- 
ground is characterized by the dimensionless quantity 



1 dp, 



GW 



p c d\nk 



(6) 



2 We have taken S„ to be 4 times larger than the contribution 
from photon shot noise alone, following Refs. [Til. |2S(I . 



III. SLOW-ROLL PARADIGM 

First, we present a Fisher matrix calculation making 
use of analytical models of slow-roll inflation. The ba- 
sic assumption here is that inflation is driven by a slow- 
rolling single scalar field and the slow-roll approximation 
is valid while the scalar field rolls down the potential from 
the point which relates the CMB scale to the point of the 
direct detection scale. 



A. Slow-roll prediction for the spectrum of the 
gravitational wave background 

In an analytic approach, it is convenient to divide the 
spectrum of the gravitational wave background into two 
parts: the initial power spectrum A 2 prim and the trans- 



fer function X/j, 



n 



GW 



1 

12 



k 
~a~H 



A? ■ 

n.prim 



{k)T 2 (k). (11) 



The initial power spectrum is predicted under the slow- 
roll approximation. In standard slow-roll inflation, a 
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scalar field <f> slowly rolls down its potential V((f), and 
the equation of motion is given as <j) + 3H<j) + V'{cf)) = 0, 
where the dot and prime denote the derivative with re- 
spect to t and 4> respectively. We define the slow-roll 
parameters in terms of V and its derivatives, 



other parameters as riT(fco) — —r/8 and ax{ko) — 
r/8 (ns + r/8 — 1). Therefore, in the framework of slow- 
roll inflation, the primordial spectrum can be written in 
terms of the parameters familiar in CMB observation, r, 
n s and A^ prim , as 



'"pi 
16tt 



>1 



'pi 



YL 

V 
V" 



8ir V 



(12) 
(13) 



where mpi = 1/vG is the Planck mass. The initial power 
spectra of scalar and tensor perturbations are predicted 
as 



Afc,, prim (&) 



L(JL 

ne \rnpi 
16 / H V 

7T \ m Pl/ 



iH 



iH 



which give the tensor-to-scalar ratio 



_ A h,prim(^) 



A? 

(, .prim 



(k) 



16e. 



(14) 
(15) 

(16) 



One may define inflationary parameters at a pivot 
wavenumber ko and express the initial power spectrum 
in a Taylor-expanded form as [IH HH 



In 



A 



ft., prim (^0 
A iprim( fc 0) 



n T (k ) In 



ko 



7-a T (fc )ln 2 T~ 
2 fco 



where the tensor spectral index nr(fc) and its running 
ar(fe) are given in terms of the slow-roll parameters 



n T (k) 

OLT(k) 



_ d^Al prim (k) 

dink 

_ driT(k) 



-2c, 



dln/c 



4erj - 8e 2 



(18) 
(19) 



Similarly, the scalar spectral index ns(k) and its running 
as(k) are given as 



n s (k) - 1 



rflnA^ prim (fc) 



rfln/c 



-6e + 277, 



a s (fc) 



_ dn s {k) 
din fc 



-16e77 + 24e 2 + 2^ 



(20) 
(21) 



where £ 2 = (mpi/27r) 4 T/'U"'/U 2 . Throughout this pa- 
per, we use parameter values evaluated at the CMB scale 
k = 0.002Mpc _1 . From Eqs. ((TBJ} and ([15]). we obtain a 
relation of single-field slow-roll inflation, which is called 
the consistency relation 



r = — 8«t- 



(22) 



From the relations of the slow-roll prediction, the 
tensor mode parameters can be connected with the 



4,prim(fc) = rA 



x exp 



In- 



£,prim 

k r 
^ + 16 



(ns + l-l) 



(23) 



The transfer function for the simple case where the 
components of the Universe are only radiation and mat- 
ter is given in Ref. jH as T 2 (k) = (3ji(fcr )/fcr ) 2 (l + 

L //„ 1„ — _ — 1 _ 

oq — 'oq — 

. The spher- 



1.34x cq + 2.5x^ q ), where x cq = k/k eq , 
6.22 x 10- 2 O m /i 2 Mpc _1 , and r = 
ical Bessel function, j\{x) = (sin a; 



x cos x) /x 2 , is re- 
placed as ji(kTo) —> l/(v2fcro) when taking the limit of 
fcro <C 1 and averaging the oscillation. Additionally, the 
amplitude of the spectrum is suppressed by the cosmo- 
logical constant [36j and changes in the effective degrees 
of freedom during the radiation-dominated era [37[ at 
direct detection frequencies. A suppression factor due 
to the cosmological constant is (1 — ^a) 2 , explained in 
the Appendix. A damping factor due to the effective 
degrees of freedom is evaluated with the temperature 
when the corresponding mode enters the horizon, Xi lc , as 
(q*(Thc)/<7*o ))(q * s o/g* s (Thc)) 4/3 , where 3*0 = 3.36 and 
g*sO = 3.90 [38[ . In the case of taking into account only 
particles in the standard model and not including SUSY 
particles or any other exotic particles, the effective de- 
grees of freedom which correspond to the direct detection 
scale are g*(Th c ) = g* s {Thc) = 106.75. Adding these two 
factors to the transfer function of Ref. [Hj], we obtain 

SI 



T 2 h {k) - (1 - n A f 



9*(Thc) 
9*0 



9*s0 

g* s {Thc) 



4/3 



V2(fcr ) 2 



(1 + 1.34x cq + 2.5a; 2 ). (24) 



We use this transfer function for the spectrum at the 
direct detection scale, which is calculated by substituting 
Eqs. (123]) and {H} into Eq. (fTTj) . 



B. Errors on inflationary parameters 

Using the analytic spectrum from the slow-roll ap- 
proximation, we calculate the Fisher matrix and fore- 
cast errors on the parameters attainable form Planck 
and BBO. We take (h, n b h 2 , n c h 2 , r, n s , r, A^ prim ) as 
model parameters. Note that the direct detection is 
not sensitive to the values of h,£lbh 2 ,£l c h 2 ,t at all. Al- 
though the transfer function given in Eq. (j24[) includes 
h, fl m h 2 (— Qbh 2 + fl c h 2 ) and 1 — Qa(= fl m ), their ef- 
fects cancel out when the transfer function is evaluated 
at higher frequencies (fc 3> k Qq ). So practically we take 
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FIG. 1: Forecast constraints on ns and r. The dotted line ellipse (red) represents the marginalized 2a confidence level by 
Planck only, and the solid one (blue) represents combined constraints from Planck and BBO. The fiducial parameters are 
'"fid = 0.1 in the left panel, = 0.03 in the right panel. 
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Variable 


Fiducial value 


% Error Planck only 


% Error Planck+BBO 


% Error CV 


% Error CV+BBO 


h 


0.724 


1.1 


1. 


0.11 


0.11 


Q c h 2 


0.108 


1.3 


1.2 


0.15 


0.15 


n b h 2 


0.227 


0.88 


0.83 


0.13 


0.13 


T 


0.089 


4.1 


4.1 


1.8 


1.8 


ns 


0.961 


0.56 


0.5 


0.11 


0.1 


r 


0.1 


29. 


16. 


1. 


1. 


A 2 

Cprim 


2.41 x 10 9 


0.79 


0.77 


0.29 


0.29 




Variable 


Fiducial value 


% Error Planck only 


% Error Planck+BBO 


% Error CV 


% Error CV+BBO 


h 


0.724 


1.1 


1.1 


0.1 


0.1 


n c h 2 


0.108 


1.3 


1.3 


0.15 


0.15 


Q b h 2 


0.227 


0.88 


0.88 


0.13 


0.13 


T 


0.089 


4.1 


4. 


1.7 


1.7 


n s 


0.961 


0.57 


0.57 


0.1 


0.1 


r 


0.03 


66. 


14. 


1.6 


1.6 


A 2 


2.41 x 10 9 


0.64 


0.64 


0.25 


0.25 



TABLE II: Marginalized lcr (68%) errors on parameters for the slow-roll inflation model. The upper table is for the rad = 0.1 
case; the lower table is for the red = 0.03 case. 



(ns, r, A? prim ) as free parameters in the calculation of 
the Fisher matrix for direct detection. Here, we show two 
cases of different fiducial values of r: r^d — 0.1 and 0.03. 
The fiducial values of the other parameters are taken to 
be the WMAP5 maximum likelihood. We use Eq. 
for tit- The gravitational wave background has a direct 
detection with SNR=18.2 in the case of r^a — 0.1, and 
SNR=8.9 in r fid = 0.03. 

Figure [T] shows the constraints in the rig — r plane ex- 
pected from Planck, and those combined with constraints 
from 10 years of observation with BBO. The Planck con- 



straints on ns and r are not particularly degenerate as 
most of the constraint for r comes from the -B-mode of 
the CMB polarization and most of the constraint for ns 
comes from the temperature measurements. The de- 
generate direction of direct detection constraints is the 
direction along which the model gives the same ampli- 
tude of the gravitational wave background spectrum at 
the direct detection frequencies. Since direct detection 
detects gravitational waves with a very narrow band- 
width (0.1 — 1Hz), it has less sensitivity to the tilt of 
the spectrum by itself and cannot measure a small de- 
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viation from the scale-invariant spectrum. Considering 
the fact that direct detection is sensitive only for the 
amplitude of the spectrum, the degeneracy line is con- 
sidered to be the direction of AQgw(/ = 0.2Hz) = 0. 
We can evaluate the degeneracy line in the ns — r plane 
assuming that A 2 prim is fixed in Eq. (|23|). which yields 
AQ GW (k) cc A(A\ pri J = [1 - rln(fc/fc )/8 + r{n s + 
r/4 - l)ln 2 (fc/fc )/16]Ar + r 2 ln 2 (fc/fc )An s /16. Sub- 
stituting fc/fco ^ 6.5 x 10 16 at k = 2ir x 0.2Hz, we 
infer that the error ellipse is elongated along the di- 
rections of 0.39Ar + 0.94An s = for r fid = 0.1, and 
0.77Ar + 0.084An s = for r fid = 0.03. This is con- 
sistent with the direction of the main axis of the error 
ellipse for Planck+BBO shown in Fig. [TJ Note that the 
error ellipses obtained by direct detection alone are much 
more elongated in these directions, but thanks to the 
tight constraint on ng from CMB, the combined ellipses 
are less elongated. Also, the direct detection constraint 
itself has no power to distinguish r and A? im both of 
which strongly affect the amplitude of the spectrum [see 
Eq. (I23p ], However, the CMB gives a quite tight con- 
straint on A 2 - prim and the degeneracy in the direction of 
A? im is broken when they are combined. 



IV. NUMERICAL CALCULATION 

In this section, we perform the Fisher matrix calcula- 
tion using the spectrum of the gravitational wave back- 
ground which is obtained by numerically solving the evo- 
lution equation of the gravitational waves. We calculate 
the evolution of the scalar field numerically and follow 
the evolution of the gravitational waves from the infla- 
tion phase up to the present. This means the amplitude 
of the spectrum obtained numerically reflects the actual 
Hubble expansion rate when each mode exits the hori- 
zon during inflation. In contrast, the slow- roll prediction 
presented in the previous section may not predict the 
precise amplitude of the spectrum at scales far from the 
CMB scale, since the spectrum is expressed by making 
use of a Taylor series approximation. However, while 
this numerical approach has the advantage of allowing 
for precise evaluation of the amplitude, this numerical 
approach requires us to assume an inflation model. Here, 
we evaluate errors on each potential model parameter of 
quadratic inflation (one potential parameter case), and 
natural inflation (two potential parameter case), which 
can give a relatively large amplitude of the gravitational 
wave background. 



As expected, we see direct detection has power mainly 
in determining the tensor-to-scalar ratio r to which CMB 
is not very sensitive. It is notable that direct detec- 
tion tightens the constraint on r more in the case of 
rfid = 0.03 as compared to rfid = 0.1. This can be ex- 
plained as follows: When ns and A 2 - prim are fixed, Eq. 
HMD yields A£1gw(^) oc [1 - r \n(k/k )/8 + r(n s + r/4 - 
l)ln 2 (fc/fc )/16]Ar. Since the direct detection does not 
have sensitivity to the frequency dependence, it is rea- 
sonable to evaluate the uncertainty fixing the frequency 
at 0.2Hz as we did above. Then the errors on r are ex- 
pected to be ay ~ Ar oc (1 - 8.5r + 23. 4r 2 )" 1 . This is 
an increasing function when r < 0.18, which means the 
error becomes smaller as r decreases. This happens be- 
cause the amplitude is not sensitive to r around r ~ 0.18 
due to the balance between the effect of r to increase 
the amplitude of the tensor-to-scalar ratio and the ef- 
fect to decrease the amplitude via the tilt of the spec- 
trum, exp(— r ln(fc/fco)/8), as seen in Eq. ([2"3"]l . There- 
fore, ^gw(/ = 0.2Hz) changes more rapidly with the 
variation of r when r is smaller, and this results in the 
smaller error on r in the rfid = 0.03 case. As r decreases 
more, oy goes to a constant value and direct detection no 
longer has power to determine r. In that case, BBO may 
give an upper limit r < 0.008 with a 3er confidence level. 



The errors on the other parameters are listed as per- 
centages of the fiducial values in Table |TTJ We see that 
direct detection does not help to improve the constraints 
on CV in this model. 



A. Method 

First, we briefly present the method of our numerical 
calculation (for details, see Ref. [Ill)- The evolution 
equation of the gravitational wave is simply expressed as 

k 2 

h* + 3Hh* + — h£ = 0. (25) 
Gr 

We evolve this equation by calculating H numerically us- 
ing the following equations. During inflation, the evolu- 
tion of the Hubble expansion is determined by the scalar 
field, which decays into radiation in the reheating phase 
following inflation. When considering a case that the de- 
cay rate T is sufficiently small, the effect of decay can be 
simply included into the scalar field equation as [4(| El[ 

4> + (3H + r)4> + V = 0, (26) 

and the energy density of the radiation p r generated from 
the scalar field obeys the equation, 

p r + AHp r = T P4> . (27) 

Then the Hubble expansion is determined by the energy 
density of the (f> field and radiation field, 

87T 

h 2 = 3^r(p* + ft0, ( 28 ) 

where the energy density of this scalar field is given as 
P<f> — < P 2 /2 + V . After the Universe become well radia- 
tion dominated, we switch to the equation for the Hubble 
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expansion rate which takes into account the change of g* 



H 2 



Hi 



g* 
9*o 



g*so 

9*s 



4/3 



n A 



(29) 

where we take Vl r h 2 = 4.15 x 10 5 and the other cosmo- 
logical parameters are set to be the WMAP maximum 
likelihood values given in Sec. Ill Al 

We derive the amplitude of the spectrum by solving 
the above equations numerically. The derivative of the 
spectrum, which is necessary for the Fisher calculation, 
is calculated by performing several calculations in which 
we change the parameters slightly. Here, not only do we 
take potential model parameters as parameters for the 
Fisher matrix, but we also take the number of e foldings 
Af as a parameter. The e-folding number is defined as 
Af(k) = ln(a e nd/afc)j where a on d is the scale factor at 
the end of inflation and is the scale factor when the 
mode exit the horizon (fc = aH) during inflation. When 
considering a mode which corresponds to the CMB scale 
ko = 0.002Mpc _1 , the value of Af is approximately given 
as 42] 



W~ 56 



2 , 10 16 GeV 
3 ln — T7T- 



1 , 10 9 GeV 
- In , 

3 Trh 



(30) 



where p cn d is the energy density at the end of inflation 
and Trh is the reheating temperature, which directly re- 
lates to the decay rate as [4(| 

THH-^ (1/4) ^) 1/4 (m P1 r)V^ (31) 



Note that we do not use Eq. (|30|) to obtain the fiducial 
value of Af , which is given numerically for a given value 
ofT. 

Here, we have assumed reheating to take place via per- 
turbative decay [H, |44| as this process can be included 
simply as in Eqs. (|26j) and fj27j> . However, in many 
cases there can be a stage of preheating (see for example 
[45l |46j]) where the inflaton decays via parametric reso- 
nance. In this paper, we choose a low decay rate consis- 
tent with decay from gravitational effects [47[ which re- 
sults in a reheating temperature of about 10 9 GeV, and 
so in the inflation models we consider there will not be a 
preheating phase. This corresponds to taking T ~ 2GeV 
in Eq. (EU). 

In order to evaluate the primordial power spectrum 
parameters, which is necessary to calculate the Fisher 
matrix of the CMB, we use the usual slow- roll formulas to 
evaluate the parameters given in Sec. 1111 Al by evaluating 
them when the mode of fc p i v exit the horizon k = aH in 
the numerical calculation of the background equations. 
This method is reasonable since the slow-roll formulas are 
sufficiently accurate for the primordial power spectrum 
at the length scales probed by the CMB. In addition to 
the model parameters taken for the CMB Fisher matrix 
in the previous section (h, tlbh 2 , £l c h 2 , r, ns, r, A^ im ), 
we also include the running of the scalar spectrum index 
as- 



B. Quadratic Inflation 

First, we investigate the case of the quadratic poten- 
tial, 



V = -mf 
2 



(32) 



We take two model parameters: the mass of the scalar 
field to and e-folds number Af, or reheating temper- 
ature log 10 (TRu/GeV) . The error in Af can be con- 
verted to the error in log 10 (TRH/GeV) by Eq. (|30|) as 
u j\f = 0.77<Ti ogio (Tk H /G G v))- The fiducial value of to is de- 
termined to satisfy the normalization of the scalar per- 
turbations A^ prim = 2.41 x 10~ 9 , and Af is determined 
by Trh = 10 9 GeV. They are derived numerically as 
(to, Af) = (1.66 x 10 13 GeV, 55.7). In this case, the grav- 
itational wave background is detected with SNR = 16.5 
by the direct detection experiment. 

Figure [5] shows the confidence contours in the m — Af 
plane, expected from Planck, BBO and both combined. 
We see Planck gives good constraints on the mass of the 
scalar field m, but has less power to determine the e-folds 
number Af. Although the BBO constraint is weaker than 
Planck, it can break the strong degeneracy in the param- 
eters since the degenerate directions are slightly different, 
and improve the errors on both of the parameters. Table 
IIIII shows how much the errors from Planck and CV de- 
crease when they are combined with the constraints from 
BBO. 

Reminding the reader that direct detection cannot dis- 
tinguish models which give the same amplitude at the 
direct detection frequencies, the degenerate direction is 
considered to be the direction in which f!cw °c A^ prim cx 
H 2 \k=aH is constant [see Eq. 1151) ]. Using the rela- 
tions that the Hubble expansion rate during inflation 
is given as H(k) 2 cx V(k) = m 2 (j>(k) 2 /2, and <p(k) re- 
lates to the e-folding number as 4>{k) 2 = 2Af(k) + 1 in 
the case of the quadratic potential model, the parame- 
ters give the same spectrum amplitude in the direction 
of m 2 [2Af(k) + 1] = const. Therefore, the degenerate 
line is considered to be Ato/to + AAf/[2Af(k) + 1] = 0. 
Substituting Af(k = 2tt x 0.2Hz) ~ 16.4 which corre- 
sponds to the direct detection scale, the degeneracy di- 
rection of the direct detection constraint is estimated as 
Am/(1.66 x 10 13 GeV) + AAA/33.7 ~ 0, which is consis- 
tent with the result shown in Fig. [5] 

We also show the constraints in terms of 
log 10 (TR,H/GeV) instead of Af. It is notable that 
the marginalized error on the reheating temperature 
is o"i og (t R h/GcV) ~ 6.0 by Planck alone, which means 
the BBN lower limit of log 10 (T RH /GeV) = -3 is less 
than 2cr away from the fiducial value without the 
direct detection constraint. When it is combined with 
direct detection constraint, the error is reduced to be 
(J iog 10 (T RH /GcV) ~ 2.8, and the BBN limit is ruled out at 
more than 4cr. 
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1.2xl0 13 1.4xl0 13 1.6X10 13 1.8X10 13 2.0X10 13 
m/GeV 

FIG. 2: Constraint on parameters of the quadratic potential. The dotted line (red) and the dashed line (green) show 
marginalized 2<r confidence region in the m — AT, log 10 Trh plane allowed by Planck and BBO, respectively, and the solid line 
(blue) shows the combined constraint. The black horizontal line shows the BBN limit on reheating temperature, Trh > lMeV. 



Variable 


Fiducial value 


% Error Planck only 


% Error Planck+BBO 


% Error CV 


% Error CV+BBO 


h 


0.724 


0.64 


0.49 


0.15 


0.15 


Q c h 2 


0.108 


0.84 


0.62 


0.2 


0.2 


Q b h 2 


0.227 


0.52 


0.5 


0.17 


0.17 


T 


0.089 


3.9 


3.7 


2.4 


2.4 


m/GeV 


1.66 x 10 13 


8.8 


4.2 


1.2 


1.1 


JV 


55.7 


8.3 


3.9 


1.1 


1.1 


log 10 (TR H /GeV) 


9.0 


67. 


31. 


9.1 


8.8 



TABLE III: Marginalized la errors on parameters for the quadratic inflation model. 




Natural inflation 



Next, we investigate natural inflation for which the 
potential is given as [48145 Oj 



V = A 



1 ± cos 



N(j) 

T 



(33) 



Here, we set 7Y = 1 and take the positive sign. In this 
case, we take three variables, A/mpi, //mpi and Af, to be 
model parameters. We investigate two cases of different 
fiducial values of /, which are taken to be / = mpi and 
/ = 2mpi, and the fiducial value of A is determined by the 
normalization of A 2 - prim . The potential becomes more 
quadratic and predicts larger amplitudes of the gravi- 
tational wave background as / increases. The fiducial 
values are, respectively, (A/mpi, //mpj, Af) = (1.07 x 
1CT 3 , 1.0, 55.2) for / = m P1 and (1.64 x 10" 3 , 2.0, 55.6) 
for / = 2mp\. The gravitational wave background is de- 
tected with SNR = 8.4 in the case of / = mpi, and SNR 
= 14.9 in the case of / = 2rapi. 

Figure [3] shows the error ellipsoids in the / = 2mpi 



case. Table IIVI gives errors on each parameters. The 
/ = 2mpi case is essentially the quadratic case and so 
the errors on Af are about the same. A Taylor expansion 
around the bottom of the potential gives the effective 
mass to be m ~ A 2 / / when / 1. Therefore the A 
and / parameters are expected to become correlated in 
the / = 2mp\ case and less so in the / = mpi case. 
This is consistent with the errors of / and A being larger 
in the / = 2mpi case. The gravitational waves are too 
low for the direct detection experiment to have much 
effect in the / = mpi case. The error in the reheating 
temperature is quite large. Even in the / = 2mpi case the 
Planck and direct detection data still has a fiducial value 
Trh = 10 9 GeV less than 3<r away from the BBN limit of 
Trh = 10~ 3 GeV. Unlike in the quadratic inflation case, 
direct detection is still useful to tighten the constraints 
from CV for the / = 2mpi model. 
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FIG. 3: Forecasted marginalized 2a constraints on parameters for the natural inflation model with / = 2mp\. CMB (red 
dotted) and CMB plus direct detection (blue solid) are shown in the f/mp\ — M and A/mpi — TV plane, respectively. In each 
case the larger contour is when the CMB is Planck and the smaller is when the CMB is cosmic variance limited. 



Variable 


Fiducial value 


% Error Planck only 


% Error Planck+BBO % Error CV % Error CV+BBO 


h 


0.724 


0.7 


0.66 


0.15 


0.15 


n c h 2 


0.108 


0.92 


0.86 


0.21 


0.21 


tt b h 2 


0.227 


0.54 


0.53 


0.18 


0.18 


T 


0.089 


3.7 


3.7 


2.4 


2.4 


A/rapi 


1.07 x 10~ 3 


8.8 


7.5 


0.57 


0.56 


f/mm 


1.0 


5.5 


4.3 


2.1 


2. 


N 


55.2 


11. 


10. 


2.3 


2.2 


log 10 (T RH /GeV) 


9.0 


85. 


83. 


18. 


18. 



Variable 


Fiducial value 


% Error Planck only 


% Error Planck+BBO % Error CV % Error CV+BBO 


h 


0.724 


0.69 


0.52 


0.15 


0.15 


tt c h 2 


0.108 


0.91 


0.66 


0.22 


0.21 


tt b h 2 


0.227 


0.53 


0.5 


0.19 


0.18 


T 


0.089 


3.9 


3.8 


2.6 


2.4 


A/rapi 


1.64 x 10~ 3 


18. 


16. 


3.3 


2.4 


//mpi 


2.0 


37. 


26. 


10. 


7.3 


N 


55.6 


9.5 


6.1 


4.1 


2.9 


log 10 (r RH /GeV) 


9.0 


77. 


49. 


33. 


23. 



TABLE IV: Marginalized la errors on parameters for the natural inflation model. The upper table is for the / = mpi case; 
the lower table is for the / = 2rapi case. 



V. CONCLUSION 



This paper is aimed at studying how direct detection of 
the inflationary gravitational wave background will deter- 
mine inflationary parameters and how it will complement 
future CMB polarization experiments. An attractive fea- 
ture of these two different methods of observation is that 
they probe two different frequencies and provide inde- 
pendent information. By calculating the Fisher matrix, 



we have investigated the degree to which the errors on 
model parameters obtained from CMB experiments will 
be reduced by direct detection of the gravitational waves 
in the BBO experiment. We have presented two different 
types of methods to calculate the Fisher matrix: One is 
evaluated analytically by making use of the slow-roll ap- 
proximation, and the other is evaluated numerically for 
the sake of more accurate predictions. In the second case, 
we also allow the temperature of reheating to be a free 
parameter. 
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In both cases, we have shown that the two different 
observations have different directions of parameter de- 
generacy and that this degeneracy is broken when they 
are combined. The degeneracy of a direct detection ex- 
periment is directed to the direction in which the model 
parameters give the same amplitude at the direct detec- 
tion frequencies. Although our result indicates that the 
BBO experiment provides a larger error in parameter es- 
timation than CMB experiments, it most certainly has 
the power to tighten constraints from Planck. We also 
found that for natural inflation, direct detection could 
even help to improve a cosmic variance-limited CMB ex- 
periment. It would be interesting to check what the im- 
provements could be in multifield inflation models where 
the number of parameters would be greater. 

Constraints on parameters which are defined by the 
slow-roll parameters were calculated by using an analytic 
spectrum of the gravitational wave background. BBO 
direct detection worked to tighten the constraint from 
Planck mainly on the tensor-to-scalar ratio. On the other 
hand, we have confirmed the intuitive fact that it does 
not have the power to constrain the other cosmological 
parameters, which makes sense since the tensor-to-scalar 
ratio is the major parameter in determining the ampli- 
tude of the gravitational wave background. Constraints 
on inflation potential parameters have been evaluated by 
using an accurate amplitude of the spectrum obtained 
from numerical calculations, for both the quadratic in- 
flation and natural inflation cases. In addition to poten- 
tial parameters, we also have taken the e-folds number, 
which corresponds to the reheating temperature, as a pa- 
rameter and have found that BBO direct detection has 
power to tighten constraints on both of the parameters. 
For quadratic inflation, we found without BBO direct 
detection, Planck could only rule out a BBN reheating 
temperature at the 2a level. However with BBO direct 
detection it could rule this out at almost the 4a level. In 
the case of natural inflation (with / = 2mpi), we found 
that a cosmic variance-limited CMB experiment could 
only rule out a BBN reheating temperature at the 4a 
level while with BBO direct detection it could be ruled 
out at almost the 6a level. 

Note that, when the signal to noise ratio is low, the 
Fisher matrix can underestimate the error bars (see for 
example [5l|). Therefore our Planck only results may 
be overly optimistic, but the combined Planck plus BBO 
and the cosmic variance CMB results should be more 
accurate. In future work, we plan to redo this analysis 
using a Markov Chain Monte Carlo method, which will 
also make it easier to add more detailed priors on the 
reheating temperature. 

For reference, let us mention another project aimed at 
detecting the inflationary gravitational wave background 
directly, called DECIGO [1(|. It has similar specifica- 
tions to BBO; the target frequency is almost the same, 
but the sensitivity is a little smaller than that of BBO. 
Because of the fact that the direction of the parameter 
degeneracy is determined only by the amplitude of the 



gravitational wave at direct the detection scale, the de- 
generate direction of constraints from DECIGO is the 
same as the one from BBO. However, since DECIGO has 
less sensitivity, the size of the error ellipse is considered to 
be bigger than BBO and the Fisher matrix analysis may 
not be applicable unless the tensor-to-scalar ratio is rela- 
tively large. (For example, while BBO detects the grav- 
itational wave background of r = 0.1 with SNR= 18.2, 
DECIGO detects this with SNR= 4.1, which may be out 
of the validity of the Fisher matrix analysis.) For this 
reason, we have presented the results only for BBO in 
this paper. 
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Appendix A: Suppression factor for the cosmological 
constant 

Reference [36) shows that the current acceleration of 
the Universe suppresses the amplitude of the inflation- 
produced gravitational waves by a factor of Q m /ri\. This 
is derived using the fact that the scale factor when the 
mode enter the horizon during the accelerating stage 
satisfies the relation <2h c oc k. However, it is not ap- 
propriate to apply this to the case of our Universe, 
f?A = 0.7, O m = 0.3, which just starts to enter the cos- 
mological constant-dominated Universe and is still not in 
the middle of a de Sitter phase of exponential expansion. 

Here, we propose a new suppression factor which gives 
a better approximation. From the behavior of the infla- 
tionary gravitational waves that keeps constant out- 
side the horizon and decreases proportional to a^ 1 after 
entering the horizon, the transfer function is considered 
to be written as Th{k) = |/ik,o|/|^k,hc| = o-Yic/o-a- The 
suppression factor is measured by comparing the trans- 
fer functions in the case of the Universe without the cos- 
mological constant, which we label with a subscript 1, 
and in the case with the cosmological constant, which we 
label with a subscript 2. Therefore, the value of inter- 
est, how much the amplitude of the gravitational waves 
is suppressed by the cosmological constant, is 

T hA k ) = fl hc,2 

Th,i(k) ahc,i 

We set the Hubble parameter Ho and the density param- 
eter of radiation £l r to be the same value in both cases. 
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The only difference is the existence of the cosmological 
constant, fiA- If we assume a fiat Universe, then the den- 
sity parameter of matter is described with the amount of 
the cosmological constant as fi m ,2 = 1 — fir — fiA, while 
fim,i = 1 — fir if there is no cosmological constant. With 
using the relation k — ah c Hi lc and rewriting the Hubble 
parameter in terms of the cosmological parameters, Eq. 
(IA1[) becomes 



Ohc,2 



k/Hi 



hc,l 



hc,l 



H J (1 — fir — fiA)a hc 3 2 + fir^hc^ + 



(A2) 



ing the matter-dominated phase. Since the contribution 
of radiation and the cosmological constant terms to the 
Hubble expansion is negligible during this phase, Eq. 
121) becomes 



(1 — fi r ) / dbcl 



a hc,l V (1 — fir — fiAj \ahc,2 



(A3) 



Neglecting the radiation density parameter, which is 
much smaller than 1, we find the suppression factor is 
approximately 



0>hc,2 



~ 1 - fi A . 



(A4) 



Let us consider a mode which enters the horizon dur- 
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